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1. Introduction 

The anomalous magnetic moment of the muon, = {g^ — 2) /2, has been both measured and 
calculated to high precision, possibly revealing a small discrepancy between Nature and the Stan- 
dard Model. The measurement of the Muon (g-2) Collaboration is a^" = 11659 208.0(6. 3)10"'° 
and has a fractional accuracy of 0.54 • 10^^. The Standard Model value has been estimated by many 
authors. One recent review |§] quotes a value of ajj" = 11659179.3(6.8)10^'°, which has just a 
slightly higher fractional error of 0.58- 10"*'. This results in a discrepancy of 3.1 a. Other theoret- 
ical estimates produce a range of discrepancies from 0.9(7 to 3.4a', but in all cases the dominant 
source of error for the Standard Model calculation is the leading order (in the QED coupling) 
hadronic contribution, a^*^. The value quoted in [^, a^'^ = 692.1 (5.6) lO"''', alone represents 
60% of the theoretical error. This quantity is a pure QCD observable and has been shown to be 
calculable in lattice QCD calculations even in Euclidean space [||]. 

This hadronic contribution to is the focus of this work. We present our initial calculation 
of a"^'^ using two-flavor maximally twisted mass fermions. This is only the second full QCD 
calculation of this quantity and represents the first such calculation to examine finite size effects 
and lattice artifacts. As we demonstrate, cleanly controlling all sources of systematic error will be 
very important to reliably calculate a^'^. 

2. Calculation 

The leading order hadronic contribution due to vacuum polarization is 



where 7^ (x) is the electromagnetic quark current. In particular, we note that the momentum integral 
in Eq. |0] is performed for > [||], thus n{c]^) can be calculated directly from lattice QCD. 
Furthermore, we remark that the momentum integral is peaked at small momentum and the kernel, 
w{q^/mj^), attains a maximal value at = (\/5 — 2)m^ 0.003 GeV^. (The inverse power of q^ 
is canceled by the subtraction n{q^) — 7r(0), which is proportional to ^^.) Meanwhile the smallest 
momentum accessible in our finite volume calculation is q^ = {In/Vf k, 0.05 GeV^. Therefore a 
reliable low q^ extrapolation is essential to calculate 7l{q^). In particular the ultra-violet subtraction 
at = required to renormalize 7t{q^) induces larger uncertainties than naively expected. 

'The range of values is determined by examining the references given in [^. The dominant source of this variation 
is the discrepancy between e+e^ and t data used to determine c^^. 




(2.1) 



where is the muon mass and w{q^/mj^) is given in [||]. The vacuum polarization, 7i{q^), is 
determined from the vacuum polarization tensor, n^viq), by 
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Table 1: Parameters used in this work. The values of a and L are given in fm and is given in MeV. 



3. Lattice Details 

We calculate n^v{^^) using dynamical maximally twisted mass fermions. The twisted quark 
mass provides an infra-red regulator that bounds the determinant of the fermion action hence elim- 
inating exceptional configurations [Q, ^. Additionally, at maximal twist physical observables are 
automatically accurate to iff{a^) in the lattice spacing [^]. 

The flavor diagonal currents retain their usual form undering twisting.^ Additionally, in the 
twisted basis we can use the conserved Noether current instead of the local current. This eliminates 
the renormalization factor required for the local current and ensures that the Ward identity holds 
even for non-zero lattice spacing. The conserved current in the twisted basis is given as follows. 



\ {z^Ta + y^.xT -xT{r- Y^)U^^X%(. } 



and has the same point-split form as the standard Wilson current. This can be understood easily 
once one realizes that both the Wilson mass term and the twisted mass term are invariant under 
local QED gauge transformations and hence do not contribute to the Noether construction of the 
conserved current. 

The calculation of n^v{q^) proceeds as for Wilson and domain wall fermions. Propagators 
from point sources at a single site and the four forward neighbors are calculated and used to con- 



struct the current-current correlator in Eq. |2.2| . The one exception is that separate u and d quark 
inversions must be performed due to the modified 75-hermiticity: YsD'lys = D^.^ 

Twisted fermions break flavor symmetry. However, the 75-hermiticity relates u and d quark 
loops and results in 7i^y{x,y) = n^*y{x,y). This expression is true for each gauge field configuration. 
The consequence is that re(7r''(^^)) = re(7r"(^^)). Hence by simply taking the real part of 7i{q^), 
which is real in the continuum limit, we eliminate any explicit flavor breaking in the valence sector. 
Additionally, we expect the real part of n to be accurate to &{a^), even if the imaginary part has 
ff{a) corrections. 

In this work we use the two-flavor dynamical gauge field configurations from the European 
Twisted Mass Collaboration [Q-^. The details of the ensembles used are given in Tab. I. Ad- 
ditionally, the hadronic contribution to has previously been calculated using quenched domain 



^This follows simply from Qjfx = exp(— i'75T30)27^ exp(— !y5T36) for 2 = 1 and T3. 

"'This can be seen from the basic loop expression y^D^^ {x,y)YvDu ' ' {y,x) = Yti^^u^ i^:y)YvY5D^^ {x,y)Y5- 
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Figure 1: Volume dependence of —n{q^). This Figure 2: Volume dependence of ^{q^) — Tt{Q). 

quantity requires an ultra-violet subtraction but is As an illustration, this quantity is renormalized at 

infra-red finite. With the exception of the lowest the lowest momentum accessible in each volume, 

point, there is no noticeable finite size effect. demonstrating the sensitive nature of the subtraction. 



wall fermions [||], quenched improved Wilson fermions [ |l0| ] and dynamical rooted asqtad improved 



staggered fermions 
4. Results 

First we examine the finite size effects in 7i{q^). We have calculated it at two volumes, L = 
2.1 fm and L = 2.7 fm, both with rriji = 310 MeV and a = 0.086 fm. In Fig. |l| we show the results 
for n{q^) for these two volumes. Although 7i{q^) is an ultra-violet divergent quantity, for a fixed a 
we expect a finite large volume limit. Fig. H demonstrates this rather clearly for all but the lowest 
point. However, the ultra-violet subtraction required to form exaggerates any differences at low 
q^. Fig. ^illustrates this point. Determining 7r(0) requires a fit, so for the purposes of illustration 
we perform the subtraction at the lowest q^ value available for each of the two volumes. Fig. ^ 
demonstrates that this subtraction can have a large effect on the resulting renormalized quantity. 
However, the integral in Eq. 2.1 is dominated by the region near q^ 0.003 GeV^ and hence is not 



fully sensitive to the overall shift in Fig. |^ 

Next we study the mj^ dependence. Figs. ^ and ^ show all five values of m^j: for a = 0.086 fm. 
In the case of wi;^ = 310 MeV, only the larger L = 2.1 fm results are shown. Fig. ^ demonstrates 
that there is no visible quark mass dependence for large q^, as expected from perturbation theory. 
Any quark mass dependence should be more visible for the low q^ region shown in Fig. 0. The 
error bars are too large to identity any quark mass dependence, however, the results do appear to 
systematically increase when proceeding from the rrijc = 450 MeV calculation down to the 270 MeV 
result. 

Now we examine the lattice artifacts in 7i{q^). We have calculated at two lattice spacings, 
a = 0.086 fm and a = 0.067 fm. In both cases we have taken m;;; = 310 MeV and L = 2.1 fm. Fig. || 
shows the unrenormalized results for n{q^) demonstrating the ultra-violet divergence present with- 
out the subtraction. In Fig. ^ we perform the subtraction, but at = 0.3 GeV^ rather than at q^ = 0. 
We see no noticeable lattice artifacts with the exception of the lowest q^ point. Unfortunately the 
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Figure 3: Quark mass dependence at large q^. There Figure 4: Quark mass dependence at low q^. There 



is no noticeable quark mass dependence at large q^, 
consistent with perturbative QCD expectations. 



is a systematic, but not statistically significant, shift 
with quark mass from m^j — 450 MeV to 270 MeV. 
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Figure 5: Lattice spacing dependence of -n{q^). Figure 6: Lattice spacing dependence of -7t{q^) 
The unrenormalized -K{q^) is shown. The discrep- matched at q^ = 0.3 GeV^. With the exception of the 
ancy illustrates the ultra-violet subtraction required lowest q^ point, there is no noticeable lattice spacing 



to renormalize n. 



dependence. 



expression for in Eq. lA requires the subtraction at = 0. This is shown in Fig. ^ where, 
as earUer, we subtract at the lowest point available in each calculation. Again, the subtraction 
induces a large difference between the results from the two lattice spacings, but we must remember 
that is dominated by values of q^ around 0.003 GeV^. 

To determine the extent to which the effects shown in Figs. ^ and ^contribute to a^^'^, we must 
parametrize and fit n{q'^) and extrapolate to = in order to perform the integral in Eq. lA. In 



Fig. |8| we show fits to polynomials in q^ with 4 terms (cubic) and 5 terms (quartic). The lattice 
results and corresponding curves are shifted vertically to illustrate the quality of the fits and the 
nature of the extrapolation to q^ = 0. Additionally, the calculation labeled m;t = 310 MeV refers to 
the larger L = 2.7 fm calculation. For all but one ensemble the cubic fit seems sufficient to describe 
the lattice results. The one exception is the niji = 310 MeV, smaller volume L = 2. 1 fm calculation 
(not shown), which requires the quartic term to accommodate the observed low q^ behavior. 
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Figure 7: Lattice spacing dependence of Ttiq^) — 
7r(0). The results have been renormahzed at the low- 
est value of at each lattice spacing to illustrate the 
effect of the subtraction. 



Figure 8: Low extrapolation. Each 7r(<7^) has 
been fit to cubic (foreground) and quartic (back- 
ground) functions of q^, showing agreement for all 
but the smaller volume at 310 MeV (not shown). 
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Figure 9: Comparison of cubic and quartic fit results 
for ajf*. There is general agreement between the cu- 
bic and quartic fits for all calculations excluding the 
smaller volume at mj^ = 310 MeV (not shown). 



Figure 10: Comparison of all full QCD calculations 
of ajf*. The results of this work are shown along 
with the staggered results of JTl] ] . Finite size effects 
and lattice artifacts are visible. 



Using the fits described above, we calculate a^^'^ using Eq. |2?T| . Fig. ^ shows the resulting 
values for the five masses at a = 0.086 fm. (Again the larger value of L = 2.7 fm is used for 
ruji = 310 MeV.) We note a clear consistency between the cubic and quartic fits as already indicated 
in Fig. p. Additionally, there is no discernible quark mass dependence as implied by Figs. |3| and 
Q In Fig. 1^, we focus specifically on the cubic fits and examine the L and a dependence of our 
results and we compare to the only other full QCD calculation |[TT|]'^. First we notice the quite 
large finite size effects and lattice artifacts as anticipated in Figs. ^ and ^ However, in general we 
find consistency with the staggered results. The agreement for the largest staggered value of rrij^ is 
quite clear. The intermediate staggered value is at a volume that is between our larger and smaller 
volumes, and the result also lies between our two results. This seems quite consistent to within the 



We will refer to the results in ||l VA as staggered and have taken the results corresponding to the cubic fits in ||l ih. 
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systematic errors that we observe near this value of mj^. Finally at the Ughtest staggered value 
of m„ appears to be a bit too high compared to our lightest value. This might be genuine quark mass 
dependence, but given the strength of finite size effects and lattice artifacts that we observe, we find 
it difficult to claim a strong quark mass dependence. In fact the sign of the discrepancy is consistent 
with both a finite size effect, which should be universal, and also with our lattice artifacts, which, 
though not universal, might still be indicative. 

5. Conclusions 

The current high precision determinations of the anomalous magnetic moment, a^, both from 
experiment and theory, indicate a small discrepancy between Nature and the Standard Model. The 
largest source of error in the theory calculation of a^j is the leading order hadronic contribution. 
We present a full QCD calculation of the vacuum polarization and, in particular, of precisely this 
hadronic contribution. We perform calculations with dynamical maximally twisted mass fermions 
with pion masses ranging from 480 MeV to 270 MeV. We observe both large finite size effects 
and lattice artifacts but find general agreement with the only other full QCD calculation. This 
work presents the first full QCD calculation of these effects and represents a first effort to begin to 
calculate controlling for all sources of error. 
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